Rapidity particle spectra in sudden hadronization of QGP 



O 

o 

(N 

C 
3 



Jean Letessier and Johann Rafelski 
Laboratoire de Physique Theorique et Hautes Energies, 
Universite Paris 7, 2 place Jussieu, F-75251 Cedex 05 

and 

Department of Physics, University of Arizona, Tucson, AZ 85721 

(June 2001) 

We show that the remaining internal longitudinal flow of colliding quarks in nuclei offers a natu- 
ral explanation for the diversity of rapidity spectral shapes observed in Pb-Pb 158A GeV nuclear 
collisions. Thus QGP sudden hadronization reaction picture is a suitable approach to explain the 
rapidity spectra of hadrons produced. 
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If a quark-gluon plasma deconfined state super cools 
and decomposes in a sudden fashion into free stream- 
ing hadrons Q, we should be able to see in the hadron 
spectra certain features of quark statistical distribu- 
tions. Such a QGP fingerprint on m^-spectra of strange 
baryons and antibaryons which are seen to be identical 
is one of key indicators for the sudden hadronization pic- 
ture §. 

A question which has remained open is how sudden 
hadronization model can be reconciled with the diverse 
shapes of rapidity spectra of observed hadrons. Great 
differences are known to be present: while hadrons made 
of quarks brought into the collision, such as, e.g., protons, 
show a pronounced double hump structure, progressively 
more peaked distributions arc seen for particles which are 
more and more composed of newly made quarks |] ^ . 

We present a schematic description of this phenomenon 
based on the assumption that the quark momentum dis- 
tribution in QGP is imprinted on hadrons produced, and 
not erased in a post hadronization period. The model we 
explore is a direct continuation of the study of m± shape 
and abundance of hadrons produced || . 

We note that many light quarks are brought into the 
collision fireball by colliding nuclei — in 158 A GeV Pb- 
Pb central interactions there are 1000+ such valance 
quark participants and there is a similar total number 
of hadrons produced. In general the two incoming quark 
fluids can retain memory of the early stage in the colli- 
sion: in the center of momentum frame two quark flu- 
ids collide with initial rapidities yi = ±2.92. In order 
to reproduce the proton data, we assume in what fol- 
lows, that the average remaining longitudinal quark flow 
is yf — ±1. This average includes quarks from qq pairs 
produced, which are mostly found at central rapidity. We 
further assume that the average flow of s,s,q quarks is 
zero. 

If quarks were not confined, their emission spectrum 
would be computed in a thermal emission model with 
flow. We proceed to establish this spectrum. Each emit- 
ting volume element is assumed to be in motion, with a 



four-vector 



and in the local frame at rest the tem- 



perature in the volume element is T, and the local quark 
distribution is of Boltzmann type. In order to count the 
states, we apply Touschek invariant phase space measure 
for the level density |?]|| : 
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Here, — Vqu^ is the flowing volume element as ob- 
served in the laboratory frame. Vo is the comoving vol- 
ume element in the local rest frame. So is the Dirac-delta 
function for the positive roots only. This is an invari- 
ant function for all proper Lorentz transformations. Us- 
ing the rapidity and transverse momentum p±_ (or mass 
m± = m 2 + variables of a particle in a cylindrical 
symmetric decomposition with respect to the collision 
axis 

p^ = (mi coshy,p± cos <f>p,p± sin<p Pl m± sinh?/) , (2) 
we recover the well known invariant measure in Eq. (hi) : 
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The four-velocity m m allowing for both transverse and 
longitudinal flows is: 



= (u ,U x ,Uy,U z ) 

uq = coshy^ cosh y || 
u y = sinh y± sin cj> u ; 

We recall the usual relations: 



u 



(4) 



= sinh y± cos <p u : 
coshyj_ sinh 



= 7-L = coshj/j_ , — = i;j_7 ± = smhyj_ , 

m m 

which illustrate the structural identity of Eq. @ with 
Eq.fl). 
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It is straightforward now to obtain the quantity re- 
quired to construct spectra Eq. (Q) , 



VP M = 7JL [ m J- cosh(?/ -y\\) - P±v± 



cos< 



(5) 



where <j> = (j) p — 4> u . For a suitable choice of coordinate 
system in which the a;-axis is pointing in the direction of 
transverse flow vector, cj> u = 0, and we can use the parti- 
cle emission angle as the azimuthal angle of integration 
for radiation emitted from the volume element: <f> = <fi p . 
This has the advantage that we can easily restrict the 
counting of particles to those emitted outwards including 
< <j> P < 7r/2 instead of the full range < 4> p < ir 

In Eq. (j^), the emitted particle variables are as usual 
y, m± and the flow vector variables used are yn and v± 
with j± = — 
explicitly is: 



The invariant spectrum written 
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There are some differences in detail between this tradi- 
tional expression and recent proposals based on con- 
sideration of the freeze-out surface geometric properties 
||. Certain features we find are the same in both ap- 
proaches: the rapidity y distribution is not appreciably 
influenced by transverse dynamics (i.e., transverse veloc- 
ity and temperature), which effectively separates. The 
m± spectra are only dependent on T, v± , but practically 
not on the longitudinal flow . We note that Eq. (j^) 
retains absolute number of particles independent of the 
flow parameters, if integral over <j) includes the full range 
< <f> < IT. 

The rapidity spectra of massless quanta in the decon- 
fined phase are shown in figure |. We see how the thermal 
spectrum (m = 0, T = 145 MeV, v± — 0.52) gradually 
'flows' apart as yii is increased from yu = (dotted) to 
yu = 1.5 (solid) in step of 0.5. Included in the double 
numerical integral (over tf> and m±) are particles emit- 
ted outwards with reference to the direction of the flow 
vector of the volume element. Model studies of kinetic 
dynamics show that inside the QGP matter we find a mix 
of fluids, the incoming quarks are still retaining some of 
the original flow while all newly made quark pairs have 
no memory of the initial condition of colliding matter. 
In particular, strange quark pairs made in the plasma do 
not flow. 

Final state hadrons are born from deconfined quarks 
with such spectra. In the previously studied case of m± 
spectra [Q it was found that the transverse flow param- 
eters of the hadronizing quark matter were in agreement 
with results of chemical particle production parameters. 
Thus the shape of m± spectra of hyperons, antihyperons 
and kaons was as if these particles were directly made 
of hadronizing QGP quanta. Addressing here the rapid- 
ity spectra we will take this point of view but proceed 



more carefully. As hyperon are formed at the fireball 
breakup, any remaining longitudinal flow present among 
fireball constituents will be imposed on the product par- 
ticle, thus A-spectra containing potentially two original 
valence quarks are stretched in y, which A-y-spectra are 
not, as they are made from newly formed particles. One 
would thus expect that anti-hyperons comprising only 
newly made quarks will appear without flow, while A- 
spectra will show flow, if such was still present. 
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FIG. 1. Rapidity spectra with flow of massless QGP 
quanta: dotted: no flow yn — 0, dashed, yn = 0.5, longdashed 
yu = 1, solid, = 1.5. 

We thus implement a particle-staged turn-on of lon- 
gitudinal flow at SPS where in central Pb-Pb collisions 
the incoming valance quark number and total number of 
hadrons produced is about the same. For protons we as- 
sume that effectively all quarks are remembering some of 
the original flow. For A with one strange quark we take 
a mix of 2/3 weight in the spectrum with flow and 1/3 
without. For particles like K + (us) we take a 50-50% mix 
and for all newly made particles like A, <I> we assume that 
only no flow components contribute. 

To describe proton rapidity spectra we assume after 
inspection of H that the longitudinal flow is y\\ ~ ±1, 
and choose this value without a further attempt to fit 
the spectra. With m q = 0, T = 145 MeV and v ± = 
0.52c (both parameters could be chosen very differently 
without altering the result, these were taken in view of 
results of Q), we obtain the proton rapidity spectrum 
we present in figure ^, bottom panel. We note that we 
have to average over positive and negative flows yu = ±1, 
since the collision in the CM frame involves both. 

The strange quark content of the A contributes with 
relative strength 33% in the central region, which results 
in a flat distribution — see the second panel from the 
bottom in figure ^. Next we take a 50-50% mix which 
gives us the result shown in the third panel from the bot- 
tom, corresponding to the K + (us) rapidity distribution 
and when resonance decays are considered, also pion dis- 
tribution. Finally, in the top panel we show prototype of 
the rapidity distribution arising for all completely newly 
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made particles such as A, <I>. Though far from perfect, 
we see in this schematic study practically exactly the be- 
havior observed. 

The same schematic rapidity behavior is seen in S-S 
collisions at 200A GeV for the particles here considered. 
A direct comparison of the charged particle distributions 
between Pb-Pb and S-S collision systems seen in Ref. 
p0{ suggests that y\\ is about 0.4 units of rapidity larger 
in the lighter collision system. Though seemingly a small 
change, this opens by 50% the gap between the fluids as 
we saw in figure [l] and a much more pronounced central 
rapidity reduction in certain particle abundances ensues. 
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FIG. 2. Rapidity particle spectra: schematic representa- 
tion within a thermal model with flow, parameters chosen for 
^/sjvjv = 17.2 GeV. See text for details. 

Our schematic method can help understand the parti- 
cle rapidity distributions at RHIC. The PHOBOS collab- 
oration has presented a nearly flat distribution of charged 
hadrons with small peaks near to pseudo rapidity r\ = 1.5 
at the edge ]iT|Jl2| . Trying our model we finds that rather 
than to have a 3-fluid quark source, we need now a con- 
tinuous distribution of sources residing in the range of 
rapidity —1.6 < y < 1.6. In other words if the hadron 
rapidity distribution is as before offering an image of the 
quark flow, we have a flow plateau as expected in the hy- 
drodynamical scaling solution proposed by Bj0rken |l3[ ]. 



In summary of this brief contribution, we have shown, 
within a schematic dynamic model, that quark flow in 
QGP can explain the diversity of particle rapidity spec- 
tra observed at ^snn = 17.2 GeV, where clearly we 
do not see yet a high energy scaling behavior. What 
we have shown is that a three component fluid, two 
made of quarks still flowing against each other, and a 
third one made of newly made components (antiquarks, 
strange quarks), is capable to imprint on the rapidity 
spectra just the right systematic behavior. This trans- 
fer of quark dynamics to hadron dynamics is similar to 
what was found in the study of the m± spectra. For 
hadrons to emerge with momentum spectra of quarks we 
need a hadronization with sudden breakup of the (super 
cooled) QGP plasma phase into free-streaming hadrons. 
Thus this study shows that such a picture is capable to 
explain rapidity spectra. 
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